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[ We calculate the cross section of Higgs boson pair production at a photon collider in the two 



Higgs doublet model. We focus on the scenario in which the lightest CP even Higgs boson (h) 
has the standard model like couplings to the gauge bosons. We take into account the one-loop 
correction to the hhh coupling as well as additional one-loop diagrams due to charged Higgs bosons 
^ \ to the 77 hh helicity amplitudes. It is found that the full cross section can be enhanced by 

' both these effects to a considerable level. We discuss the impact of these corrections on the hhh 

coupling measurement at the photon collider. 
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The Higgs sector is the last unknown part of the standard model (SM) for elementary par- 
ticles. Discovery of the Higgs boson and the measurement of its properties at current and 
future experiments are crucial to establish our basic picture for spontaneous electroweak 
symmetry breaking (EWSB) and the mechanism of particle mass generation. The Higgs 
mechanism would be experimentally tested after the discovery of a new scalar particle by 
measuring its mass and the coupling to the weak gauge bosons. The mass generation mech- 
anism for quarks and charged leptons via the Yukawa interaction is also clarified by the 
precise determination of both the fermion masses and the Yukawa coupling constants. If 
the deviation from the SM relation between the mass and the coupling is found, it can be 
regarded as an evidence of new physics beyond the SM. The nature of EWSB can be revealed 
through the experimental reconstruction of the Higgs potential, for which the measurement 
of the Higgs self-coupling is essential|l|, [2, ^, 4, ^. The structure of the Higgs potential 
depends on the scenario of new physics beyond the SM^, Q], so that the experimental de- 
termination of the triple Higgs boson coupling can be a probe of each new physics scenario. 
Furthermore, the property of the Higgs potential would be directly related to the aspect 
of the electroweak phase transition in the early Universe, which could have impact on the 
problem of the electroweak baryogenesisQ]. 

It is known that the measurement of the triple Higgs boson coupling is rather challeng- 
ing at the CERN Large Hadron Collider (LHC), requiring huge luminosity. A study has 
shown that at the SLHC with the luminosity of 3000 fb^^, expected accuracy would be 
about 20-30% for the mass (m^) of the Higgs boson (h) to be around 170 Gev|l|, Q]- At 
the international linear collider (ILC), the main processes for the hhh measurement are 
the double Higgs boson production mechanisms via the Higgs-stlahlung and the W-boson 
fusion 3|, IJ]. If the collider energy is lower than 1 TeV, the double Higgs strahlung process 
e~^e~ — ^ Zhh is important for a light Higgs boson with the mass of 120-140 GeV, while 
for higher energies the W-boson fusion process e^e~ hhvv becomes dominant due to 

n 

its t-channel nature [5|. Sensitivity to the hhh coupling in these processes becomes rapidly 
worse for greater Higgs boson masses. In particular, for the intermediate mass range (140 
GeV < rrth < 200 GeV), it has not yet been known how accurately the hhh coupling can be 
measured by the electron-positron collision. 

The photon collider is an option of the ILC. The possibility of measuring the hhh coupling 
via the process of 77 — > hh has been discussed in Ref. ^|, where the cross section has been 
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calculated at the one-loop level, and the dependence on the triple Higgs boson coupling 



constant is studied. In Ref. 10| the statistical sensitivity to the hhh coupling constant 
has been studied especially for a light Higgs boson mass in relatively low energy collisions. 
Recently, these analyses have been extended for wider regions of the Higgs boson masses 
and the collider energies. It has been found that when the collision energy is limited to be 
lower than 500-600 GeV the statistical sensitivity to the hhh coupling can be better for the 
process in the 77 collision than that in the electron-positron collision for the Higgs boson 
with the mass of 160 GeV 111 ]. 



Unlike the double Higgs production processes e~^e~ Zhh and e~^e~ — * hhuu in e^e~ 
collisions, 77 — hh is an one-loop induced process. When the origin of the shift in the 
hhh coupling would be due to one-loop corrections by new particles, it may also affect the 
amplitude of 77 hh directly through additional one-particle-irreducible (IPI) one-loop 
diagrams of 77/1 and '-ffhh vertices. 

In this letter, we consider the new particle effect on the 77 — ^ hh cross sections in the 
two Higgs doublet model (THDM), in which additional CP-even, CP-odd and charged Higgs 
bosons appear. It is known that a non-decoupling one-loop effect due to these extra Higgs 
bosons can enhance the hhh coupling constant by (9(100) % 6|. In the 77 — > hh helicity 
amplitudes, there are additional one-loop diagrams by the charged Higgs boson loop to the 
ordinary SM diagrams (the W-boson loop and the top quark loop). We find that both the 
charged Higgs boson loop contribution to the 77 — > hh amplitudes and the non-decoupling 
effect on the hhh coupling can enhance the cross section from its SM value significantly. 
We consider how the new contribution to the cross section of 77 — > hh would affect the 
measurement of the triple Higgs boson coupling at a 77 collider. 

In order to examine the new physics effect on 77 — > hh, we calculate the helicity ampli- 
tudes in the THDM. We impose a discrete symmetry to the model to avoid flavor changing 
neutral current in a natural way 12| . The Higgs potential is then given by 

^THDM = /ii|$l|^ + /i2l$2p- (/i3'^l$2 + h.C.) 

+Ai|$i|^ + A2|$2|^ + Aal^iH^sr + A4|$I$2r + y {{^{^2? + h.c.} , (1) 

where $1 and $2 are two Higgs doublets with hypercharge -1-1/2. We here include the soft 
breaking term for the discrete symmetry by the parameter In general, yU.| and A5 are 
complex, but we here take them to be real for simplicity. We parameterize the doublet fields 



as 



00, 



IZi) 



1,2), 



(2) 



(246 GeV)2. The 



where vacuum expectation values ( VEVs) Vi and f 2 satisfy vl + v"^ = v 
mass matrices can be diagonahzed by introducing the mixing angles a and /5, where a diag- 
onalizes the mass matrix of the CP-even neutral bosons, and tan/5 = V2/V1. Consequently, 
we have two CP even {h and if), a CP-odd {A) and a pair of charged {H^) bosons. We 
define a such that h is the SM-like Higgs boson when sin(/9 — a) = 1. We do not specify the 
type of Yukawa interactions 3], because it does not much affect the following discussions. 

Throughout this letter, we concentrate on the case with so called the SM-like limit [sin(/5— 
a) = 1], where the lightest Higgs boson h has the same tree-level couplings as the SM Higgs 
boson, and the other bosons do not couple to gauge bosons and behave just as extra scalar 
bosons. In this limit, the masses of the Higgs bosons are^ 



H 



{Ai cos^ (3 + \2 sin^ /5 + 2(A3 + A4 + A5) cos^ (3 sin^ (3}v'^ 
M^ + \ {Ai + A2 - 2(A3 + A4 + A5)} (1 - cos4/3)?;2, 



2 



(3) 
(4) 
(5) 
(6) 



where M(= |/i3|/-\/sin P cos P) represents the soft breaking scale for the discrete symmetry, 
and determines the decoupling property of the extra Higgs bosons. When M ~ 0, the extra 
Higgs bosons H, A and receive their masses from the VEV, so that the masses are 
proportional to A^. Large masses cause significant non-decoupling effect in the radiative 
correction to the hhh coupling. On the other hand, when M ^ v the masses are determined 
by M. In this case, the quantum effect decouples for M — > 00. 

There are several important constraints on the THDM parameters from the data. The 
LEP direct search results give the lower bounds rrih > 114 GeV in the SM-like limit and Trifj, 
m^, fTijj± ~ 80-90 GeV 14]. In addition, the rho parameter data at the LEP requires the 



approximate custodial symmetry in the Higgs potential. This implies that ?77,^± ~ or 
sin(/9 — a) ^ 1 and mH± — ttlj^. The Higgs potential is also constrained from the tree level 



^ For the case without the SM-Hke Hmit, see Ref. [7| for example 
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FIG. 1: Feynman diagrams for charged Higgs boson contributions to 77 hh. Wavy lines 
represent photons, dotted hnes in loops represent charged Higgs bosons H^, and other dotted lines 
are the neutral Higgs bosons. 



unitarity jisl. . the triviality and vacuum stability 17|, in particular for the case where the 
non-decoupling effect is important as in the discussion here. For M ~ 0, masses of the extra 
Higgs bosons H, A and are bounded from above by about 500 GeV for tan/3 = 1, when 
they are non-zero M, these bounds are relaxed depending on the 

value of M. The constraint from b s'-f gives a lower bound on the mass of if^ depending 



on the type of Yukawa interaction; i.e., in Model II ?Ti^± > 295 GeV (95% CL) 18 1 
Recent data for B ru can also give a constraint on the charged Higgs mass especially for 
large values of tan/3 in Model Iljigl. [2^. In the following analysis, we do not include these 
constraints from B-physics because we do not specify type of Yukawa interactions. 

In the THDM with sin(/3 — a) = 1, the one-loop helicity amplitudes for the initial photon 
heUcities ii and £2 (^i = +1 or — 1) are given as 

where Xhhh = ~3m^/t;, A4{ii, (^2, ^hhh) is the SM amplitude given in Ref. j^, and 
AA4{ii, £2, Xhhh) represents additional one-loop contributions from the charged Higgs boson 
loop to the 77 hh cross section. We note that Xh^h has the same form as in the SM when 
sin(/3 — a) = 1. Due to the parity we have A1thdm(-^15 -^2) = -^thdm(— ^1, — ^2), so that 
there are independent two helicity amplitudes. 

The Feynman diagrams which contribute to AA^ are shown in Fig. 1. AM. is given for 
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each helicity set for sin(/3 — a) ^ 1 as 

s [4 J 

-1 {in - ml) {Dir + + + } , (8) 

and 

a^AM{+^ -, A...) = -{\H.H-r^^ {in - ml) {Dlf^ + + + ) , (9) 
where s, t and u are ordinary Mandelstam variables for the sub processes, and 



ijkl 
ab 



), C'o(s) = Coiml,ml, s, mj^±, mj^±, m^±), and D 
Dghipf, p'j^p t^ Ph m^±,m^±,m^±,m^±). Here we employ the Passarino-Veltman formalism 
in Ref. 21 1. We take the same normalization for these amplitudes as in Ref. j^. We note 
that AA^(+, — , Xhhh) is independent of \hhh because of no s-channel diagram contribution. 
The scalar coupling constants and X^hH+H- defined by 

The relative sign between A1(£i,£25 ^hhh) and AA1(£i,£25 ^hhh) has been checked to be con- 



sistent with the results for the effective Lagrangian in Eq. (19) in Ref. 22| in the large mass 
limit for inner particles. 

In Eq. (I7j), Xhhh is the tree level coupling constant. It is known that in the THDM 
Xhhh can be changed by the one-loop contribution of extra Higgs bosons due to the non- 
decoupling effect (when M ~ 0). In the following analysis, we include such an effect on the 
cross sections replacing Xhhh by the effective coupling r™^°^(s, m^, m^), which is evaluated 
at the one- loop level asQ| 



pTHDM/^ 2 2^ 

hhh '"'h^ '"'h) — 



ml / Nrm 



19-0-21,2^2 \ 



<S>=H,A,H+,H- 



\2'K'^v'^m'i V m| J Sir'^v'^ml 



As a striking feature, there are quartic power contributions of the masses of extra Higgs 
bosons which is divided by v'^ml, when M ~ 0. Thus, the large mass of the extra Higgs 
boson {H, A, H^) with the lighter SM-like Higgs boson h would cause large quantum 
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corrections to the hhh coupling, which amount to 50-100%. This effect can be regarded as 
the leading two loop contribution to 77 — > hh in our analysis. The exact one-loop formula 
for r™^^^ is given in Ref. [3], which has been used in our actual numerical analysis. 
Finally, the cross section for the each subprocess is given by^ 

-^r~ ~ ^I-^™dm(^i,4)| , (12) 
where A^thdm(^1' ^2) is defined by 

-Mthd^(^i,^2) = -M(^i,4,r™) + AM{i^,i,,Tlj^n■ (13) 

We comment on the consistency of our perturbation calculation. One might think that 
the inclusion of the one-loop corrected hhh vertex function T™^^^ in the calculation of the 
cross section 77 hh would be inconsistent unless we also take all the other two loop 
contributions into account. Our calculation can be justified in the following sense. First 
of all, r™^^^ is a gauge invariant subset. Second, it can be seen from Eq. (ITT]) that the 
deviation from the SM value Ar™DM/rSM ( r™DM/rSM _ ^^ere F^m is one-loop 
vertex function of hhh in the SM given in Ref. [7|, can be of 0{1) for the case of M^, m\ <C 
m|,, whereas the contributions from the other two loop diagrams do not contain the factor 
m|,/m^, and thus relatively unimportant for m| ^ m^. Therefore, we can safely neglect 
these effects as compared to the non-decouphng loop effect in the hhh coupling. The details 
are shown in Appendix. 

In Fig. [21 the cross sections of 77 — hh for the helicity set (+, +) are shown as a 
function of the photon-photon collision energy E^^. In the left [right] figure, parameters 
are chosen to be nih = 120 GeV [nih = 160 GeV], sin(/3 — a) = 1, tan/3 = 1, M = and 



400 GeV. In this case, AF™°M/rSM amounts to about 120% for 
ruh = 120 GeV (80% for nih = 160 GeV)[6| ^. The five curves in each figure correspond to 
the following cases, 

(a) THDM 2-loop: the cross section in the THDM with additional one-loop corrections 
to the hhh vertex, F)^™^; i.e., the contribution from A1thdm(+; +) ^'i- p3l) . 



^ The right hand side of Eq. ([T2l) is different from the formula in Ref. Q by factor 1/2, but Eq. ([121 



reproduces figures shown in Ref. 
3 The results of ATlJ^^^yrl^,^ with A4 ^ are given in Ref. 
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FIG. 2: The cross section &{+,+) for the sub process 77 — > hh with the photon hehcity set (+, +) 
as a function of the colhsion energy E^^. In the left [right] figure the parameters are taken to be 
ruh = 120 [160] GeV for m^(= rrijj = = m^^) = 400 GeV, sin(/3 — a) = 1, tan/3 = 1 and 
M = 0. 



(b) THDM 1-loop: the cross section in the THDM with the tree level hhh coupling con- 
stant Xhhh] i-e., the contribution from A^thdm(+' +) ^'i- ©• 

(c) SM 2-loop: the cross section in the SM with additional top loop correction to the hhh 

n 

coupling r^^^ given in Ref. [7|. 

(d) SM 1-loop: the cross section in the SM with the tree level hhh coupling constant A^J^ 
(= Xhhh for sin(/3 -a) = 1). 



(e) For comparison, we also show the result which corresponds to the SM 1-loop result 

^THD 
hhh 



with the effective hhh coupling F™^^ 



In the left figure, there are three peaks in the 2-loop THDM cross section. The one at the 
lowest E^^ is the peak just above the threshold of hh production. There the cross section 
is by about factor three enhanced as compared to the SM prediction due to the effect of 
AF™^/FSM 120%) because of the dominance of the pole dia grams in 77 hh. The 
second peak at around E^y^ ~ 400 GeV comes from the top quark loop contribution which 
is enhanced by the threshold of top pair production. Around this point, the 2-loop THDM 
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cross section in the case (a) can be well described by that in the case (e). For E^^ ~ 400-600 
GeV, the cross section in the THDM 2-loop result deviates from the SM value due to both 
the charged Higgs loop effect in and the effect of AFj^'^^/F^^^. The third peak at 
around E^^ ~ 850 GeV is the threshold enhancement of the charged Higgs boson loop in 
AA^, where the real production of charged Higgs bosons occurs. The contribution from the 
non-pole one- loop diagrams is dominant. In the right figure, we can see two peaks around 
E^^ ~ 350-400 GeV and 850 GeV. At the first peak, the contribution from the pole diagrams 
is dominant so that the cross section is largely enhanced by the effect of AVj^^^ /T^f^ by 
several times 100% for E^^ ~ 350 GeV. It also amounts to about 80% for E^^ ~ 400 GeV. 
For E^^ < 600-700 GeV, the result in the case (e) gives a good description of that in the 
case (a). The second peak is due to the threshold effect of the real H^H~ production as in 
the left figure. 

The full cross section of e~e~ —^'yy^hh is given from the sub cross sections by 
convoluting the photon luminosity spectrum[9|]: 



da = / dr—^ —^da{+, +) + -^^da{+, -) , (14) 

where ^/s is the centre-of-mass energy of the e~e~ system, and 

dL^^ _ Z"^" dy 



—f^{x,y)f^{x,T/y), (15) 

Jr/yr^ y 

where t = s/s, y = Ey/Eb with E^ and Ef, being the energy of photon and electron beams 
respectively, and y^ = x/(l+x) with x = AEbUOo/m?^ where ujq is the laser photon energy and 
me is the electron mass. In our study, we set x = 4.8. The photon momentum distribution 
function fy{x,y) and mean helicities of the two photon beams {i = 1,2) are given in 



Ref. [23 1 . 



In Fig. [3l the full cross sections of e~e~ ^ '-ff ^ hh are shown for rrih = 120 GeV in 
the left figure and = 160 GeV in the right figure, respectively, as a function of ^/s for 
various values of the extra Higgs boson masses (= rrijj = uij^ = m^j^) in the cases of 
tan/? = 1, sin(/5 — a) = 1 and M = 0. In order to extract the contribution from a{+, +) 
that is sensitive to the hhh vertex, we take the polarizations of the initial laser beam to be 
both —1, and those for the initial electrons to be both +0.45 [9|. The full cross section for 
= 400 GeV has similar energy dependences to the sub cross section (t(+, +) in Fig. [2l 
where corresponding energies are rescaled approximately by around y/s ~ -E'^^/0.8 due to 
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FIG. 3: The full cross section of e e 



77 



hh as a function of y/s for each value of m^(= 



H 



m , 



m 



with sin(/3 - a) = 1, tan/3 = 1 and M = 0. The case for ruh = 120 [160] GeV 



is shown in the left [right] figure. 



the photon luminosity spectrum. For smaller m^, the peak around i/i ~ 350 GeV becomes 
lower because of smaller AF™°^/F^Jf;,. 

In Fig. m the full cross sections are shown as a function of for rrih = 120 GeV at 
y/s = 350 GeV (the left figure) and rrih = 160 GeV at ^/s = 600 GeV (the right figure). In 
each figure, five curves correspond to the cases (a) to (e) in Fig. 2. The other parameters are 
taken to be sm{(3 — a) = 1, tan f3 = 1 and M = 0. In the left figure, one can see that the cross 
section is enhanced due to the enlarged F™^^^ for larger values of which is proportional 
to m| (when M ~ 0). This implies that the cross section for these parameters is essentially 
determined by the pole diagram contributions. The effect of the charged Higgs boson loop 
from AM is relatively small since the threshold of charged Higgs boson production is far. 
Therefore, the deviation in the cross section from the SM value is smaller for relatively small 
(10-20% for m$ < 300 GeV due to the charged Higgs loop effect in AA^) but it becomes 
rapidly enhanced for greater values of (O(IOO) % for m$ > 350 GeV due to the large 
AF™^^^). A similar enhancement for the large m$ values can be seen in the right figure. 
The enhancement in the cross section in the THDM can also be seen for < 250 GeV, 
where the threshold effect of the charged Higgs boson loop in AAi appears around ^/s ~ 600 
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FIG. 4: In the left [right] figure, the full cross section of e e — > 77 ^ /i/i at \/s = 350 GeV 
[600 GeV] for rrih = 120 [160] GeV is shown as a function of m^{= rrijj = = m^±) with 
sin(/3 - a) = 1, tan/3 = 1 and M = 0. 



GeV in addition to that of the top quark loop diagrams in Ai. For = 250-400 GeV, 
both contributions from the charged Higgs boson loop contribution and the effective hhh 
coupling are important and enhance the cross section from its SM value by 40-50%. 

We have analysed the new physics loop effects on the cross section of 77 hh in the 
THDM including the next to leading effect due to the extra Higgs boson loop diagram in 
the hhh vertex. Our analysis shows that the cross section can be largely changed from the 
SM prediction by the two kinds of contributions; i.e., additonal contribution by the charged 
Higgs boson loop in AA4, and the effective one-loop hhh vertex T'^j^j^^ enhanced by the non- 
decoupling effect of extra Higgs bosons. The cross section strongly depends on rrih and ^/s 
and also on m^. The approximation of the full cross section in the case (a) (2-loop THDM) 
by using the result in the case (e) (SM-|-F™^°^) is a good description for a/s <^ 2m$/0.8. 
On the other hand, in a wide region between threshold of top pair production and that of 
charged Higgs boson pair production, both the contributions (those from AAi and from 
^h™^) important. In the region below the threshold of the real production of extra 
Higgs bosons, the cross section can be a few times 0.1 fb in the THDM while that in the SM 
is about 0.05 fb. Such differences from the SM prediction would be detectable at a future 
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photon collider. 



24| in which extra 



We note that the analysis in this letter can be applied to the models 
charged scalar bosons appear with a potentially large loop correction in the hhh coupling. 
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Note added: After this work was finished, we noticed the paper [25|] which studied 77 — >■ 
hh in the THDM. Our paper includes the additional contribution of the hhh vertex (the 
leading two- loop effect on 77 — > hh), which was not considered in {25! . 



Appendix 

If the mass of the particle in the loop comes from the VEV, a large mass implies a large 
coupling constant, so that a naive argument of the decoupling theorem is not applied. It is 
known that in such powerlike mass contribution of particles in the loop appears in 

the one-loop contribution. This is called the non-decoupling effect. 

When one-loop corrected hhh vertex Tj]^^^ largely deviates from F^^^ due to the non- 
decoupling property of the extra Higgs bosons, the main two loop contribution to 77 — > hh 
comes from the s-channel diagrams with the effective hhh coupling. We here show this by 
the use of a power counting method. For simplicity, we consider the leading powerlike effect 
of the mass of particles in the loops in the two-loop diagrams for the case with M ~ where 
masses of extra Higgs bosons are proportional to the VEV so that the non- decoupling effect 
is maximal. 

When M ~ 0, the coupling constants of hH^H~ and hhH^H~ are proportional to 
m^±/t> and m^±/t;^, respectively. We consider the situation that rnjj± ^ y/s > 2mh. The 
leading non- decoupling effect of the one-loop triangle-type diagram in Fig. [H^ up-right) 
and that of the one-loop box-type diagram in Fig. [Tt^bottom-left) are evaluated as 



167r2 V s — m\ \ v J (Airv) 
1 

r\j 

167r^ ^2 (47rw)" 
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FIG. 5: Example of the two-loop diagrams contributing to 77 — > hh. 

where we used the fact that the effective 77/1 and '-f-fhh vertices come from the dimension 
six operator so that they are proportial to q'^/v and g^/f^ at the leading order, 

respectively, where s. Therefore, the effect of mH± on 77 — * hh can be at most 

logm^± at the one- loop level. A similar conclusion of power counting can also be obtained 
for one-loop effects of top and bottom quarks and W bosons to 77 — ^ hh. 

Next, let us examine two- loop diagrams shown in Fig. [51 The non-decoupling effect in 
the diagram (a) in Fig. [5]( up-left) is calculated as 



M 



2— loop 
(a) 
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(18) 



where momenta of external lines are neglected, and k is the momentum in the loop of the 
effective hhh vertex, which is replaced by the greatest dimensionful parameter of the system; 
i.e. This result of the power counting is not changed even after the renormalization 

of the hhh vertex is performed 6|. There are other two loop diagrams which are generated 
from the s-channel type one-loop diagram, such as the diagram (b) in Fig. [S]^ up-right) where 
there is the bridge of h in the triangle type loop. Its non-decoupling effect is evaluated 
as 
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.(19) 



The dependence on fn^^± is not quartic but quadratic. We have examined all the other 
two-loop diagrams which are generated from the one-loop s-channel diagram and confirmed 
that they are the same or less power dependence on as the diagram (b). 

A similar counting can also be applied for the diagrams such as the diagram (c) in 
Fig. down-left) where charged Higgs bosons are running in the both loops, and the diagram 
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(d) in Fig. [5]( down-right) where ladder of h is added to the one-loop box type diagram; 



We find that all the IPl two-loop diagrams of ^^hh also have the quadratic or less power 
dependences on 'mii±. 

The power dependence on mj^^i in the two point function of h can be reduced by the 
renormalization of mass m^, but the highest power of rtiii± in the IPI two loop diagrams of 
'-yjhh does not change by the renormalization. 

In conclusion, the non-decoupling effect of on the renormalized amplitude of 77 hh 
at the two loop level can be described as 

^ 77^2 + ^ (jT^) + ^ (tT^) 1 ' (22) 

where the second term in RHS comes from the s-channel diagrams which include the one-loop 
corrected hhh vertex. For the case where non-decoupling property of the extra Higgs bosons 
is important, the contribution from this term is dominant when ttLj^^ ^ m/,. Although we 
gave the explanation for the charged Higgs loop effects, this argument can also be applied 
to loop effects of all quarks, gauge bosons and extra Higgs bosons with non-decoupling 
property. 
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